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Abstract. We consider previously well-known models in epidemiology 
where the parameter for incubation period is used as one of the impor- 
tant components to explain the dynamics of the variables. Such models 
are extended here to explain the dynamics with respect to a given therapy 
that prolongs the incubation period. A deconvolution method is demon- 
strated for estimation of parameters in the situations when no-therapy and 
multiple therapies are given to the infected population. The models and 
deconvolution method are extended in order to study the impact of ther- 
apy in age-structured populations. A generalisation for a situation when 
n— types of therapies are available is given. 

Key words: Epidemic models, deconvolution, conditional probability. 
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1. Preliminaries, basic model and modified models 

The incubation period is generally defined as 'the time duration between the 
time a virus or bacteria enters the human body and the time at which clinical 
clinical symptoms occur'. This duration can vary from case to case depending 
upon the route through which the virus or bacteria enters the immune system 
of an individual and in some cases depends upon the age of the infected indi- 
vidual. For influenza this duration is 1 - 2 days, for common cold 2-5 days, 
for measles 8-12 days, for SARS a maximum of up to 10 days, for rubella 
14 - 21 days, and for HIV infection to AIDS 6 months to 10 years or more. 
The incubation period can be used as a measure of rapidity of the illness after 
interaction with the virus or bacteria. It is not easy to collect information on 
the incubation period of infected individuals unless they are monitored. One of 
the direct ways of estimating the average incubation period of a given virus in 
the population is by surveillance and followup of the infected individuals until 
they develop symptoms. All the infected individuals may not be aware of their 
infection until symptoms appear and followup is subject to the availability of 
an individual. It might not be possible to follow up individuals in a typical 
situation, where time taken for the onset of symptoms from the infection is 
longer, especially if infected individuals are lost to followup. Hence, there are 
limitations on directly estimating the everage incubation period from prospec- 
tive cohort studies. Nevertheless, the incubation period occupies an important 
role along with other parameters in modelling the disease spread and under- 
standing the basic reproductive rate. A useful description of various epidemic 
models, and of estimation of parameters like the incubation period, transmis- 
sion rates, forces of infections are presented in Anderson and May (1991). The 
degree of importance of obtaining accurate average incubation periods varies 
with the incubation period of the disease. This degree of variation causes 
mathematical models to act sensitively in predicting future burden. Models 
describing dynamics of disease spread where the incubation period is shorter 
are less subject for producing misleading results than models for the spread 
with longer and varying incubation periods. Especially for predicting AIDS, 
the epidemic models developed, depend heavily on parameters that determine 
transmission rates of infection from infected to susceptible and on the pa- 
rameter which explains the average time to progress to AIDS. A review of 
various modeling approaches and quantitative techniques to estimate the in- 
cubation period can be found in Castillo-Chavez (1989) and Brookmeyer and 
Gail (1994). A set of parametric and non-parametric models are proposed to 
estimate the incubation period when the data is doubly-censored (Rao et al 
2005). The introduction of anti-retroviral therapies and protease inhibitors 
during the 1990s in several parts of the world resulted decline in opportunistic 
infections related to AIDS (Mouton et al (1997) et al (1999), Conti, S et al 
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(2000), Hung et al (2003)). As a result of such interventions, the incubation 
period was prolonged. There have been attempts to estimate the incubation 
period that vary due to drug intervention using statistical density functions 
(Artzrouni (1994)). The impact of this variation on the HIV dynamics, stabil- 
ity and on basic reproduction number has been investigated (Castillo-Chavez 
et al (1989)). In this section, we first consider a simple ODE model that ex- 
plains the dynamics of HIV spread in a population leading to AIDS. We then 
consider a similar model where incubation period is a variable with respect 
to a given therapy. We address issues of estimating incubation period to be 
used in such dynamical models and the impact of above mentioned therapies. 
Various ideas and the outline of this work are given at the end of this section. 

Perhaps the most fundamental model for the epidemiology of AIDS is that 
given by Anderson and May (1988, 1991), which takes the form 

(1.1) ^ = XX-(d + fi )Y, 

Here the total population (N) is divided into susceptibles (X), infectives 
(Y) and individuals with the full blown disease (D z ). The parameter A is the 
input into the susceptible class, which can be defined as the number of births 
in the population, A is the force of infection, \i is general (non-AIDS related) 
mortality, 7 is disease related mortality and 1/d is the average incubation 
period. Here the incubation period is defined as the duration of time between 
infection and onset of full blown disease. 

There has been an observed increase in the mean incubation period since the 
availability of therapies for AIDS (UNAIDS, 2002). Drugs are available which 
cannot eliminate virus from the body, but are helpful in prolonging the life 
of an individual by slowing the disease progression (in other words increasing 
the incubation period). For example, Protease inhibitors (say drug 1) facil- 
itate in producing non-infectious virus (only infectious virus participates in 
new virus production), hence slowing the disease progression; anti-retroviral 
therapy (say drug 2) blocks virus from interracting with the non-infected cells 
and hence reduces the infection process within the cell population (see section 
5, Nowak and May (2000) and Perelson and Nelson (1999) for fuller details); 
and a combination of the above two drugs (say drug 3) can be more effective 
by simultaneously combining the function of drug 1 and drug 2. Note that, 
when model (1.1) was developed, the above described drugs were not available. 
We assume that once individuals start taking drugs, their average incubation 
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period is prolonged. So, instead of assuming a constant 1/d, we assume that 
it varies based on the drug type. Thus we define 1/di = J R Zig(zi)dzi, for 
i = 0,1,2,3, where i = denotes the without drug scenario, % — 1 for drug 
1, i — 2 for drug 2 and i = 3 for drug 3. Here g is the probability density 
function with a certain parameter set (say B) and Z{ is a continuous random 
variable representing the incubation period. Here Zi is a real valued function 
defined on a standard probability space (S, A, P), where S is the space of 
elementary events, A is called a Borel fields, and P(A) is probability of the 
event A e A. So, z : S — > R. We can also denote this integral as a Stieltjes 
integral LzidG(zi), where G{z) = P(Z < z). We further assume without loss 
of generality that 



[1.2) / z dG(z ) < / Zl dG{z x ) < z 2 dG(z 2 ) < / z 3 dG(z 3 ), 

/ z dG(z ) < / Z\dG(zx) > / z 2 dG(z 2 ) < / ^3^(2:3). 
Jm. Jr Jr Jr 



(In the next section, we will give a detailed estimation procedure for B.) 
Applying these varying incubation periods, model f 1 1 - 1 h is modified as follows: 
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dt 
dD 



dt 



X X- UJ^z dG(z )J +V>Y , 
XtX-lf [ z^GizS] +i\y x , 



dX 

— = A-(A + Ai + A 2 + A 3 + /i)X, 
dt 

dYv 

dt 

dYx 
~dt 

dt 

dYs 
dt 

dP Z0 
dt 




X 2 X-U / z 2 dG(z 2 ) +v}Y, 



-1 



X 3 X-{[ / z 3 dG(z 3 ) +v}Y 3 , 



-i 



/ z dG(z ) I Y - j D ZQ , 
Jm / 



-i 



;i-3) ^ = ( / 



/ ^(2^ Y 3 - l3 D z:i , 
Jm. / 



where Y ,Y 1 ,Y 2 and F 3 are variables for infectives, D ZQ , D Zl , D Z2 and D Z3 
are variables for individuals with the full blown disease, A ,Ai,A 2 and A 3 and 
7oj 7ij 72 and 73 are variables for disease related mortality for no-drug, drug 1, 
drug 2 and drug 3 respectively. General mortality and disease-related mor- 
tality are incorporated in to the model to demonstrate the basic structure of 
the model, and our aim here is to estimate B and thus to estimate JL ZidG(zi) 
for all i such that simulations of the model are performed. In model (ll.3p . the 
total population N = X + Y +Y 1 +Y 2 + Y 3 + D Zo + D Zl + D Z2 + D Z3 satisfies 



i=0 i=o 

Estimation of parameters for the varying incubation periods is important 
for understanding the impact of drugs in prolonging the onset of disease and 
thus to prolong the life. The set B will also be useful in obtaining varying 
basic reproductive rates, for all i = 0,1,2,3. This can be computed as 
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Roi = M» / K ZidG(zi) by assuming independence of the impact of various drugs. 
So far, there is no evidence that f3, the probability of infecting a susceptible 
partner, changes with the activation of a drug in the body. If we assume this 
as a constant, then i? o > -Roi > or < R Q2 > Ro3- Roi < {Roi, R02, R03} , 
because individuals are assumed to have longer incubation period due to the 
effect of drugs. In the absence of clinical evidence, we assume that the impact 
of drugl and drug2 follows any one of the following relations: J R zidG(zi) < 
or > f R Z2dG{z'i). Similarly, another important epidemiological measure, the 
doubling time, tdi is obtained as t<a = hi(2) J R ZidG(zi)/ [Roi — 1] . 

Anti-retroviral therapy helps in blocking the virus from interacting with 
cells and simultaneously providing protease inhibitors helps in producing non- 
infectious virus. So without loss of generality, it is assumed that the impact of 
double drug therapy is better than a single drug therapy. If we assume disease 
related mortality is constant for all % then 7$ = 7. Where there are n types of 
drugs available, we write the general form for the above dynamical model as 
follows: 



As a special case we can consider all the parameters in the above model 
as Stieltjes integrals and can estimate them using the rigorous procedure ex- 
plained in the next section. 





ni 



(1.4) 
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This paper is organised as follows: In section 2, we describe in detail the 
estimation of the set B for up to three drugs; section 3 gives the corresponding 
expressions for the conditional probabilities of N— drugs. We construct some 
theoretical examples in section 4 to demonstrate the method explained in 
section 2. There were several attempts to study, through ODE models, the 
role of anti-retroviral therapy and protease inhibitors in controlling the virus 
multiplication inside the human body. In section 5 we briefly review these 
models. In section 6 analysis for age-structured populations is described in 
detail. Overall conclusions are given in section 7. Appendix I gives equations 
for conditinal probability when incubation period for various drug types does 
not have the monotonicity property. Appendix II gives some more theoretical 
examples when the incubation period is truncated to the right. 

2. Varying incubation periods and conditional 



In this section, we will give a detailed procedure to estimate B through a 
deconvolution technique. Let B be split into a collection of four parameter 
sets say, B={_B , Bi, B 2 , B 3 } for the four types of scenarios described in the 
previous section. Let H be the time of infection and Z be the incubation 
period, then the time of onset of the disease can be represented as D = H + Z. 
There have been studies (see for list of references Brookmeyer and Gail (1994)), 
in which H and Z were assumed independent and D was estimated through 
convolution. We outline the general idea of convolution and then give the 
convolution of H and Z. Suppose (a n ) and (b n ) are two sequences of numbers 
over the time period, then 



where (a n ) * (b n ) is the convolution of these sequences with an operator '*'. 
Suppose a and b are mutually independent random variables and let Ac(x) 
and Bc(x) be their Laplace transformations, then a + b has the Laplace trans- 
formation AcBc- Since the multiplication of the Laplace transformation is 
associative and commutative, it follows that (a n ) * (b n ) is also associative and 
commutative. Instead of discrete notation, suppose a and b are continuous 
and independent with probability density functions h and g, then the density 
of h + g is given by 



PROBABILITIES 



n 



(2.1) 




k=0 



/CO POO 
h(t - s)g(t)dt = / h(t)g(t-s)ds 
■CO J —CO 



RUNNING HEAD: VARYING INCUBATION MODEL 



9 



Suppose G(s) = g(s)ds, and F(s) = f(s)ds then 



(2.2) F(s) = / h(t)G(t - s)ds. 

J —oo 

We call F the convolution of h and G. Suppose the above h and G represent 
the infection density and incubation period distribution function; then the 
convolution of h and G represents the cumulative number of disease cases 
reported (or observed), and is given by 



/oo 
h(t)G(t - s)ds. 
■oo 



The kind of convolution in (12.31) was used to estimate the number of AIDS 
cases for the first time by Brookmeyer and Gail (1988). Information on G may 
not be available for some populations. In such situations, G has been estimated 
through deconvolution from the information available on h * G and h ( Rao 
and Kakehashi, 2004, 2005). In this section we will construct conditional 
probabilities for each drug type and express the function that maximizes B. 
These kind of conditional probabilities derived for the drug type were not 
available earlier for the incubation periods when the total number of reported 
disease cases were considered. Note that h * G is the cumulative number of 
disease cases. 

Let X , Xi,X 2 , X n _k, X n _i, X n __ m , ...,X n be the disease cases avail- 
able in the time intervals [£4-1, E/j) for i — 0, 1, 2, ...n — k, ...n — /, ...(n — m) + 
1, ...n+ 1. Suppose E is the event of diagnosis of disease after the first infection 
at T . Let E = {E , Ex, E 2 , E 3 } and E occurs in the interval [Z7 , Z7„_^) 5 E\ 
(or E 2 ) in [C/ n _ fc , U n _i) (or [C/ n _j, U n - m )) and E 3 in [U n ^ m , U n ). Now D, the 
cumulative number of disease cases up to time U n , can be expressed from (12.31) 
as follows: 



D(U <s<U n ) = / h(t)G(t - s)ds + / h(t)G(t - s)ds 

Jo Ju n _ k 



rU n - m pU n 

(2.4) + / h{t)G{t - s)ds + I h(t)G(t-s)ds 

•J Un — l J U n — m 



h(t)G(t - s)ds + \ 

'U n -l JU n 

D (A, B/U n ) = jf n - k h{t/A )G(t - s/B )ds + h{t/A 1 )G{t - s/B x )ds 
+ h(t/A 2 )G(t - s/B 2 )ds + f£_ m h(t/A 3 )G(t - s/B 3 )ds. 
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In the above equation A ,Ai, A 2 and A3 are the parameter sets for the h for 
drugO, drugl, drugl and drug3. An infected individual could fall in to any 
of the intervals described above, and similarly a full-blown disease diagnosed 
individual could fall in the same interval, but for a given individual the chrono- 
logical time of infection would be earlier than that of diagnosis of the disease. 
U n -k is the time of introduction of drugs after infection at U$. Individuals 
who were diagnosed on or after U n - k , and before U n , were taking one of the 
three drugs. If E x G [U n - k , U n -i) and E 2 G [U n _i, C/ n _ m ) Z 1 < Z 2 , otherwise 
if E l G \U n - h U n -m) and E 2 G [U n - k , U n - t ) then and if E u E 2 G [U n - k , C/ n _ m ) 
then Z\ = Z 2 . An individual who was diagnosed with the disease before 
U n must have developed symptoms in one of the four intervals [Uo,U n - k ), 
[U n - k ,U n -i), [U n -i,U n - m ) and [U n - m ,U n ). Let Ej G [Uj-^Uj) C [U ,U n - k ), 
then the conditional probability of the occurence of Ej given E is expressed 
as 



P (Ej/E) = P (Uj-! <D< Uj/D < U n ) 

D{A Q ,B Q /U j )-D{A Q ,B Q lU j . 1 ) 



D(A ,B /U n ) 



h(t/A )G(t - s/B )ds. 



h(t/A )G(t - s/B )ds 



(2.5) 



u. 



2-1 



h(t/Ao)G(t - s/B )ds. 



h(t/A )G(t - s/B )ds 



If drugs were initiated at then these conditional probabilities con- 

structed above will change according to the occurence of Ei, E 2 ,E 3 . Consider 
E 1 n E 2 = 0. Let E k G [U k -i, U k ) C [U n _ k , [/"„_,), and E x G [U n _ k , [/"„_,), then 

P(E k /E) = P{U k ^<D<U k /D<U n ) 

D(A 1 ,B 1 /U k )-D(A 1 ,B 1 /U k _ 1 ) 
D(A 1 ,B 1 /U n ) 



U k 



hit/AJGit - s/B x )ds. 



h(t/A x )G(t - s/BJds 



fUk-i r rU„ 

(2.6) -/ h(t/A 1 )G(t- s/B 1 )ds. / h(t/Ai)G(t - s/B-^d. 
Jo Uo 



-1 
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Suppose E k G [U k -i,U k ) C [U n - k ,U n -i), and E 2 G [U n - k , U n -i) i.e a situa- 
tion when Z\ > Z 2 , then 



P{E k /E) = / h(t/A 2 )G(t-s/B 2 )ds. 
Jo 



/ h(t/A 2 )G(t - s/B 2 )ds 
Jo 



(2.7) 



rUk-i 
Jo 



h(t/A 2 )G(t - s/B 2 )ds. 



h(t/A 2 )G(t - s/B 2 )ds 



Let Ei G [Ui-i,Ui) C [C/ n _ z , C/„_ m ), and £ 2 G [[/„_;, C/ n _ m ), then 



P = P (U1-1 < D < Ui/D < U n ) 



D{A 2 ,B 2 /U l )-D{A 2 ,B 2 /U l _ 1 ) 



D{A 2l B 2 /U n ) 



h(t/A 2 )G(t - s/B 2 )ds. 



h{t/A 2 )G{t - s/B 2 )ds 



(2.8) 



Ui- 



h(t/A 2 )G(t - s/B 2 )ds. 



U„ 



h{t/A 2 )G{t - s/B 2 )ds 



IJo 



Suppose Ei G [Ui-i,Ui) C [U n -i,U n - m ), and Pi G [C/ n _; , U n - m ) i.e a situa- 
tion when Z\ > Z 2 , then 



P(Ei/E) 



u. 



hit/AJGit-s/BJds. 



-, -1 



h(t/Ai)G(t - s/B^ds 



(2.9) -/ h(t/A 1 )G(t- s/B 1 )ds. / h(t/A 1 )G(t- s/Bjd 

Jo Uo 

Now consider E 1 = E 2 G [C/ p _i,C/ p ) C [[/„_&, U n - m ), i.e. Zi = Z 2 , then the 
conditional probabilities contain the same parameter sets. In this situation, 



P (E p /E) 



hit/AJGit-s/BJds. 



u„ 



h{t/A 1 )G{t - s/B 1 )ds 



(2.10) -/ h{t/A 1 )G{t-s/B 1 )ds. \ h{t j A x )G{t - s j B x )dt 

Jo YJo 

Since Z 3 > Z Q , Z x , Z 2 , suppose E 3 G [U m -i, U m ) C [C/ n _ m , C/ n ], then 
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P{U m .!<D <U m /D <U n ) 



D (A 3 , B 3 /U m ) - D (A 3 , B 3 /U m ^) 



D(A 3 ,B 3 /U n ) 
Therefore, 



P(E m /E) 
(2.11) 



u„ 



h{t/A 3 )G{t - s/B 3 )ds. 



u n 



h{t/A 3 )G{t - s/B 3 )ds 



—i 



h(t/A 3 )G(t - s/B 3 )ds. 



Un 



-, -1 



h{t/A 3 )G{t - s/B 3 )ds 



The above conditional probabilities P (Ej/E) , P (E k /E) , P (E t /E) , P (E p /E) 
and P (E m /~E) are the probabilities associated with the intervals [£/f- 1, Uy) 
, [Z7fc/_i, Up), [U V ^U V ), [Upi^Upr) and [Z7 m /_i, ?7 m ') for the ranges of j,k, I, p 
and m defined above. Since, Xo, Xi, X 2 , X n _fc, X n _;, X n _ m , X n are 
mutually exclusive, we assume they follow a parametric distribution with the 
above probabilitiesare mutually exclusive, so we assume they follow the multi- 
nomial property of the distribution of the values in the time intervals and the 
above conditional probabilities. Then the likelihood functions corresponding 
to the event set E are L (A,B/P,-) =nj=i Pj (A,B/2», £1(2) (A,B/P fc ) = 
K~=n-k p k' (A, B/2V), L m (A, B/P,) = rirln-^ fl' (A, B/l», L 1=2 (A, B/P p 



nn-m 
v'=n- 



P p , (A, B/7» and L 3 (A, B/P„ 



1 lr 



„ P ra * (A, B/T„ 



Here 



P. = P(P./E). We estimate A by fitting an infection curve from the inci- 
dence data and we then estimate B by maximizing the likelihood functions 
expressed above. The best estimate of A could be information for initial 
values of X and Y in the model (11.3 1) . Using the corresponding estimate 
of B, we obtain J R ZidF(zi). In such situations, the above likelihood func- 



tions would be L 
Li 



nr" p T] 



2 = njjju p p pl and L 3 



- n n 

1 lm'=n-i 



n— Z 



P 



re; 



P 

n—k k 



-2(1) 



re: 



:n— Z Z 



3. Generalization for multiple drug impact 

In this section, expressions for the conditional probabilities are presented 
when multiple drugs are administrated in the population. Let N = {N , Ni, N 2 , N N } 
be the number of available drugs and Z = {Z , Z\, Z 2 , Z^} be their corre- 
sponding incubation periods. Further let Z < Z\ < Z 2 <,...,< Z^ and A 
and B be their parametric sets. Then 
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rU N() 

D(A,B/U Nn ) = / h(t/A No )G(t-s/B No )ds 
Jo 

+ [ Nl h(t/A Nl )G(t - s/B Nl )ds 

J Un 

I-Un n 

(3.1) ■■• + / h(t/A NN )G(t - s/B NN )ds. 



Now, P{E N JE) = P {U Ni _, <D< U N JD < U Nn ) and L Ni (for some i) 
can be computed as follows: 



P(EnJE) 



h(t/A N .)G(t - s/B Ni )ds x 



U Nn n - I 

h(t/A NN )G(t-s/B NN )ds 



(3.2) 



/ h{t/A Ni _AG{t - s/B Ni _ x )ds x 
h{t/A NN )G{t-s/B NN )ds 



o 







rij^jVi-j J ^ s maximized for the set [Aj, Sj] by the procedure ex- 
plained in the previous section. We will obtain N sets of [A, B] values, and 
the corresponding likelihood values are -La^, L^ 2} L^ 3 , L^ N . In the above, 
we have assumed monotinicity of [Zi] to arrive at ( 13.2I ). If the [Zi] values are 
not monotonic then the various conditional probabilities can be constructed 
as explained in the previous section. There we explained the general ex- 
pression when there were a finite number of drugs available on the market. 
A detailed construction of various conditional probabilities is not necessary 
for the purpose of the present section (for details see Appendix I). When 
the ZiS are not monotonic, and if they follow some order, say for example, 
Z Q > Z\ < Z 2 > ... < Z N , then the conditional probabilities can be con- 
structed in the same way as equations (12.7H2.9I) were. Suppose (Z p ) are equal 
for each p, then there will be two scenarios arising: one for before drug in- 
tervention and one after drug intervention. For this situation, the likelihood 
equation is L Np = n^Ar p _i Pp P where P ( K E Np / E) , is given as follows: 
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P(E N JE) = / h(t/A Np )G(t-s/B Np )d. 



S X 



Jo 



1 



h(t/A NN )G(t-s/B NN )ds 




h(t/A Np _AG{t - s/B Np _Ads x 



(3.3) 



- rU N 

Jo 



N 



h(t/A N „)G(t- s/B NN )ds 



4. Theoretical examples 



In this section, we show some examples of the likelihood function constructed 
in the previous section, to estimate A, and B. Let h(s) follow a quadratic 
exponential and B follow a) a gamma function, and b) a logistic function. 
Infections in most of the countries started declining after the availability of 
antiretroviral therapies (see Conti, S et al (2000), Hung, C-C et al (2003)), and 
incidence in the recent period was found to be stable in some countries like 
India (Rao et al (2006)). This motivated us to choose a quadratic exponential 
to represent h(s), namely h(s) = exp(ais 2 +a2S+as) for all — oo < a\, a 2 , 03 < 
00. A quadratic exponential function has been shown to be a good model for 
representing the above declines in the incidence rates (see Rao and Kakehashi, 
2005). The incubation period for AIDS is large as well as variable, therefore, 
functions like the gamma, Weibull and logistic can mimic several shapes to fit 
the incubation period data depending on their parameter values. Such well- 
known functions were used by many researchers for modeling the incubation 
period of AIDS. We now demonstrate the application of such functions for the 
theory explained in section 2. 

4.1. Example 1: Gamma function. If u > is the parameter and T(u) 
is the complete distribution function, then the incomplete gamma distribu- 
tion, G(uj; tj) = p^jj e~ x x UJ ~ 1 dx, for a > 0, tj > and a + tj ^ 0. From the 
conditional probability equations from (12.50 to (12.111) . and the likelihood equa- 
tions explaind in the later part of section 3, the following are the likelihood 
equations without a drug and for with three types of drugs: 



(4.1) L (a 1 ,a 2 ,a 3 ;uj/P j ) = ai(j)a 2 (j) - JJ a x {j - l)a 2 (j) 

j j 
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where 



ai(J ~ 1) 



Too 



Uj — s) w duj > ds 



u 3-l 



ais +02S+03 



1 f j 



tj-i 



x <j — / e - ^"- 1 — ^ttj-i - s) w-1 d«j_i [> (is 



a 2 (j) = 



yo 



eais 2 +a2S+Q3 ^ l_ 1^ e _ (un _ s) ^ _ s y-i dUn (, ds 



(4.2) Li( 2 ) (a 1 ,a 2 ,a 3 ;uj/P k ) = Y[ a 1 (k)a 2 (k) - a x (k - l)a 2 (k) 



where 



a x (k) = 
a^k-1) = 



u k 



.,01^+023+03 



ruk-i 

.1 



X 



]_ ftk-l 



e- (u fc-i- 3) (M fc „i - s) w ~MM fc _i } (is 



a 2 (k) 



OlS 2 +02S+03 



rw 



e -(.u n -s)^ Un _ s y L du n >ds 



(4.3) L 2 (i)(a 1 ,a 2 ,a 3 ;uj/Pi) = j J ai{l)a 2 (l) - ai(Z - l)a 2 (l) 



where 
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di(Z) 
ai (l-l) 



v-l-i 



a 2 (l) = 



Jo 



ais 2 +a2S+ci3 



X < ft So e ~ (n ~ 1 ~ s) ( u '-i - S ) W 1<iu '-i rfs 



(4.4) L 1=2 (a 1 ,a 2 ,a 3 ;uj/P p ) = JJai(p)a 2 (p) - JJai(p - l)a 2 (p) 



where 



ai(p) 
ai(p- 1) 



eQas 2 +a2S+a3 \1_ f P e -iu P -s) ( Up _ s Y'^du p \ ds 







3 ais^+a2S+a3 



02(P) = 



/•tin 

yo 



x <j — J e-0-i- s »(Mp-i - s)^ 1 ^! [> ds 



(4.5) L 3 (ai,a 2 ,a 3 ; w/P m ) = \\ a 1 (m)a 2 (m) - ai(m - l)a 2 (m) 



where 
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ai(m) 
ai(m — 1) 



a 2 (m) 



-I 



t m \ -\T„ 

e -(u m -s) ( Mm _ s) UJ ~ 1 du m I ds 



«m-l 



^ rt m -i 1 i Trr 



1 /"*" 

rw 



4.2. Example 2: Logistic function. Suppose 9i,9 2 are parameters and 
F(9 1 ,9 2 ;t j ) = < 1 + e _(J ^ _) I , for 9 1 ,9 2 > 0, is the distribution function. 



The likelihood equations to obtain the parameters of logistic distribution with- 
out drugs and for three types of drugs are as follows: 



(4.6) L (a 1 ,a 2 ,a 3 ;9 1 ,9 2 /P j ) = JJ a[(j)a 2 (j) - JJ a[(j <- l)a' 2 (j) 

3 3 

where 



a' 2 (j) = 



I e aiS +a2S+c *3 1 J 1 + 9 2 ' 



i ii i u j-T-~ e t\ 

D ais z +a2S+a3 ) ) ^ _|_ ^ ) 



,a 1 s 2 +a 2 s+a; i J J ^ _|_ g -( x ) 



Ti 



ds 



Jo 



ds 



(4.7) L 1(2) (a 1 ,a 2 ,a 3 ;9i,9 2 /P k ) = J~\ a'i{k)a' 2 (k) - JJ a\(k - l)a' 2 (k) 



where 
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a[(k) = 

a\{k-l) = 
a' 2 (k) = 



a 1 s 2 +a2S+a 3 ) ) ^ _|_ g _ ( — e~ 2 — ) 



i n 



ds 



Jo 



a 1 s 2 +a 2 s+a 3 ) J ^ _|_ g -( k g 2 — -) 



ds 



<H*-+<>:>*+<>:! J J ]_ _|_ e ~( U "e 2 t ' 1 ) | 



ds 



(4.8) L 2(1) (ai,a 2 ,a 3 ;0i,02/fl) = l\^(l)a' 2 (l) -- JJa'^l - l)a' 2 (l) 



where 



a[(l) = 

a[(l-l) = 
a' 2 (l) = 



Ul 



a 1 s 2 +a 2 s+a 3 ) ) j _|_ g -( g 2 ) 



(is 



cns 2 +a2S+a 3 ) J ^ _|_ e _ ( _ ^ _ ^ — ~) 



(4.9) L 1=2 (ai,a 2 ,tt3;^i,^2/ J Pp) = f\a'i(p)a' 2 (p) -\\_a\iyp - l)a' 2 (p) 



where 



a[(p) = 

a'i(p-l) = 
a' 2 (p) = 



ais 2 +a 2 s+t«3 J J ^ _|_ g _ ( P 9 2 ) 



ds 



g ais 2 +a 2 s+a3 ^ | ]_ _|_ g^C'e/ 1 ) 



ds 
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(4:.10j, 3 (a 1 ,a2,a 3 ;9 1 ,9 2 /Pm) = ai(m)a' 2 (m) - a'^m - l)a 2 (m) 



where 



a[ (m) 
a'-yim — 1) 
a' 2 (m) 



I | / "m-l- 9 l \ 



1 + e v «2 



-1 



ds 



(is 



4.3. Example 3: Log-normal function. Suppose and cr are parameters 

and LNF(fi,a;tj) = \ jl + er/ (^f) } , for //, <r > 0, is the distribution 

function. (Here er/{.} is the error function of the Gaussian function). The 
likelihood equations to obtain the parameters of the logistic distribution with- 
out drugs and for three types of drugs are as follows: 

(4.11) L (ai,a 2 ,a 3 ;/W^) = II a "UKU) " II <^ ~ ^KU) 



where 



°>"U) 



1 f Uj 



= 







,ais 2 +a 2 s+a a J 1 _|_ er/ 







In a; — /i 
\nx — /j, 



a 



V2 



ds 



ds 



ds 



(4.12) L m a 2 , a 3 ; ^, <r/P fc ) = ]J a'l(k)a%(k) - ]J a'/ (A; - l)a%{k) 



where 
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a'[{k) 
a'l(k - 1) 
a»(k) 



-u fe 



eai sl+a 2 s+a 3 ) i + er f 



v>k-i 



e c n s 2 +a 2 s+a 3 J ^ _|_ er j 



lnx — fi 
ax/2 



In a; — ji 



ds 



eai sl+a 2 s+a 3 } 1 + er f 



ax/2 



lax — fi 
ax/2 



-\T k 



ds 



ds 



(4.13) L 2{1) {a u a 2 ,az-n,<r/Pi) = \{ a'[ (l)a' 2 ' (I) - J[ a'[{l - l)a' 2 ' (I) 

i i 

where 



<(/) 
a'[{l-l) 
a' 2 \l) 



; ais2+a2S+ft3 \l + erf 



o 

1 r^- 1 

2 Jo 

1 r 

2 J 



lux — fi 
ax/2 



ds 



eaiS i +a2S+a3 \ \ ■ < rf i hl^lL ] = (/ , 



ax/2 



Ti 



(4.14) L 1=2 ( ai , a 2 , a 3 ; /i, <r/P p ) = J] <(pK(p) - J] <(p - l)a' 2 '(p) 

p p 

where 



<(P) 
<(p-l) 

a' 2 '(p) 



1 /-Mp 
27o 







ax/2 



1 

2J 



e a lS * +a2 s + a 3 ) 1 + erf { lnX J 1 ] K (h 



ax/2 



(4.15)L 3 («!, a 2 , a 3 ; cr/P m ) = Q ai(ra)a' 2 (m) - a'^ra - l)a 2 (ra) 

m m 

where 
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a'[{m) = 
a"(m-l) = 



5. Models for Drug efficacy 

In the past 10-12 years there have been models developed to study HIV 
dynamics in vivo and the role of drugs in changing the virus dynamics inside 
the human body. In this section we outline such models given by others and 
explain how specific drugs intervene the standard dynamics (For details on the 
models given in this section see Ho et al (1995), Wei et al (1995), Perelson and 
Nelson (1999), Ferguson et al (1999) and Nowak and May (2000)). 

Let Xi, Fiand Z 1 represents variables corresponding to uninfected T cells, 
infected T cells and HIV virus in the serum. Then the system explaining the 
dynamics of these variables is given by 



a' 2 '(m) = 



^ = M-^X-il-^PiXZ 
dY 

(5.1) = {l-U^XZ-^Y 

^ = M^Y x -c x Z x 
at 

Here the parameter for reverse transcriptase drug (i.e. anti-retroviral ther- 
apy) is £r G [0, 1]. A value of indicates no effect of the drug in blocking the 
virus from infecting an uninfected cell, and 1 indicates that the drug is block- 
ing the virus with 100 per cent effectiveness. Ai is the rate of generation of new 
target cells, \i x is the death rate of these uninfected target cells and fi\ is the 
rate at which virus particles infect these uninfected target cells. Once the cells 
are infected it is assumed that they will die at a rate \i\ due to the immune sys- 
tem or the virus. This is based on the assumption that the probability of cell 
death at time t is given by an exponential distribution with an average lifetime 
of 1/ 1 \i x for uninfected cells and for the infected cells. It is also assumed 
that an infected cell will produce M new virus particles during their lifetime, 
so this number multiplied by the average life time gives the average rate of 
virus production i.e. M\i\. Lastly, c\ is the constant clearance rate at which 
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virus particles are cleared from the system. The steady state solutions for the 
model (jSHJ) are X* = Ai//^, Y* = 0, Z\ = and X{ = ci/M/3i(l-£ K ), Y* = 
(A x //ii) (1 - 1/Ro), Z{ = (A 1 M / 3 1 (1 - U) - V-xCl) /Pi(1-Zr)*- The basic re- 
productive rate for the model (15.11) is i?o = \\M(3\ (1 — £r) j \i x c\ and will be 
nullified when reverse transcriptase is successful. The solution will be stable 
when R > 1. 

When protease inhibitors are introduced we will have four variables in the 
system. Infected cell will produce infectious virus and non-infectious virus due 
to drug intervention. Let Zj and Z NI be the variables corresponding to infected 
and non-infected virus. If fa and (3% are the infection rates, then (3iX\ Zj is the 
mass action term and (3% = due to drugs. Before giving the drug, whatever 
infectious virus is there will be there for some time and infect other cells. In 
the blood plasma there will be a concentration of infectious virus particles Zj 
and noninfectious virus particles Zjvj (Z = Zi + Zm). Quantitatively, assume 
that 100 percent effective PI is given, then it will not have any impact on 
producing infected T cells from the virus, but it will have impact on infected 
T cells and block production of Zj, and only Z^i is produced. And whatever 
infectious virus remaining in the blood plasma will be there for a short time 
and will be cleared exponentially, i.e. Zj(t) = Z e~ Clt . Then the basic model 
for the virus dynamics becomes 



(5.2) 



dXi 

~df 
dYx 

~dT 
dZi 

~dt 
dZ NI 

dt 



Ai - /i x Xi - (3 2 X x Zi - o{f3zXxZ N i) 

fcX 2 Zj - 

-c x Z 1 



M/nYt - Cl Z 



Nl 



The assumption of 100 percent efficacy in protease inhibitors is not realistic, 
so an effectiveness factor, £p is introduced in the system (Perelson and Nelson, 
1999) as follows 
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Ai — \i x X\ — P2X1Z1 



(5.3) 



dt 
dYi 



dt 
dZi 



(1 - 60AO/iKi - Cl z z 




— ciZ NI 



dt 



The steady state solutions for the model (15.31) are X* = \\/fi x , Y* = 



0, = 0, X* = Cl /M/3 2 (l-ep), >? =(AiM) (1 - l/i?o), Z} = (Ai/ Cl ) 
(1 — ^ P ) M — (^//3 2 ) and = AiM£ P / Cl (1 - 1/iZo). Here i? = AiM/3 2 



(1 — £p) I \i x c\. The basic reproductive rate is mathematically not different 
for RT and PI. In other words, the reduction in R will be same as that of 
efficacy due to RT intervention. Herz et al (1996) discuss the alternatives 
of prescribing both anti-retroviral therapy and protease inhibitors simultane- 
ously. We assumed this double therapy is more effective than single theray in 
our analysis presented in sections 2 and 3. 



In this section we extend the models 1X731 and fTT4l to accommodate age struc- 
ture into the population mixing and epidemiology parameters. The incubation 
period for children is shorter than that of adults. Within the adult popula- 
tion there could be variability due to age at the time of infection. There are 
studies that analyse the HIV data on age collected at the time of infection 
to study parameters like incubation period (Becker et al, 2003), and some 
studies incorporate age structure in the models to explain the impact of an 
age-dependent incubation period (Griffits, 2000). In the absence of availability 
of cohort data, the methods explained in section 2 could be of great use to es- 
timate the incubation period. The analysis and method explained there could 
be done for individuals of each age (single ages). However, we describe the 
age-structure model and the method to obtain the incubation period in this 
section by considering j age groups. In a hospital set-up it is relatively easy 
to follow cohorts of age groups compared to following cohorts of individuals 
for each age group. 

Suppose the population in the j th age group is divided into Xj susceptible, 
Y Q j, Yij, Y 2t j, Y 3 j are infected and D ZQi j, D Zl j, D Z2i j, D Z3 j individuals with 
the disease without drugs, and for drugl, drug2, drugS respectively. The 
differential equations explaining these variables are 



6. Age-structured populations 
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(6.1) 



dX 3 

dt 
dt 

dt 

dt 

dY 3J 
dt 



— Aj - (A° fe + \] k + \ 2 jk + X^ k + fj,j) Xj, 



dD 



20 ,3 



dt 
dD 



zi,3 



dt 
dD 



22 J 



dt 
dD 



23 J 



dt 



\° jk Xj 



^)k X 3 



tfk X 3 



tfk X 3 



z jdG(z j) + fij > Y 0jj , 



zijdG(z ltj )) +Hj>Y 1J , 



z 2 , j dG{z 2 j)\ +iij>Y 2J , 



z 3 ,jdG(z 3:j )) +fi j }Y 3J , 



z 0J dG(z 0d ) Y 0J - 7 0ii D. 



20 ,3 



zijdGizu) Y ld - lhj D 



21 J 



dF(z 2tj )\ Y 2 j- l2d D 



/ Z %3 

Jr 

/ z 3tj dG{z 3d ) 
Jr 



-l 



z 2 ,3 



Ys, 3 ~ TfsjD. 



3^Z3,3- 



Here, Aj is the input of susceptibles for the individuals in the age group 
j, fj,j is the mortality rate, A° fc ,A] fc ,\ 2 k and A? fe are the forces of infection at 
which a susceptible in the age group j is infected by an infected individual in 
the age group k and7oj, 71 j, 72 ,j and 73^ are disease related mortality rates 
for the infected individuals in the age group j without drugs, and with drugl, 
drug2, and drugS for the individuals. {f R Zi,jdG(zij)') is the rate of disease 
progression for the infected individual for the age group j for the drug type i. 

If there are n drug types available then the general model describing the 
dynamics of various variables described above is as follows: 
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dt 
dt 



/ z o,jdG(z ,j) Y j - 7o,jD Zod , 
Jr J 



(6.2) 



dD z n ,j 

dt 



\ z n jdG(z n j) | Y n> j "fn,jD Zn j 
Jr J 



where a^j is the mortality rate of infected individuals of drug type % in the 
age group j. 

6.1. Varying incubation periods for age-structured populations. We 

are interested in the average incubation period for a group of individuals in the 
age group j. If H(j) is the time of infection and Z(j) is the incubation period 
for the j th age group, then the time of onset of the disease for this age group is 
D(J) = H(j) + Z(j). This is the time of onset of the disease for an individual 
who acquired the infection while in the j th age group. Development of the 
disease will be some time units (for example: months, years) after infection at 
age j. An individual who acquired the infection at age j is assumed to develop 
the full disease before completion of the same age j or > j. Given H(j), 
for some j, then D(j) is allowed to occur at age j'(f — j, j + 1, • • • ,j + cu, 
where j + u is the last age group for the possibility of infection). Clearly, 
H{j) < D(J). H(j) = D(j) is possible if an individual acquired infection 
and attains disease before completion of age j. One can do analysis using a 
bi-annual (or half-yearly) ageing process. 

Consider an infection and disease development matrix (see figure 2) where 
each cell (j,f) denotes the (infection age groups, disease onset age groups) for 
j = 0, 1, 2, j + iv ; j' = 0, 1, 2, j + iv. Only those cells for which j < j' are 
provided, and other cells are left blank for which the incubation period is not 
defined. In the matrix, all the eligible cells are denoted, so obviously there are 
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D(j) 

1 ... j j+1 ... j+w 



j+w 



(0,0) 


(0,1) 




(0,j) 


(0,j+l 




(0,j+v 




(1,1) 




(l,j) 


(1 j+1 




(l,j+w 






















(jj) 


Oj+i) 




(j,j+w) 










a+ij+i) 




(j+lj+W i 




























(j+w,j+w) 



Figure 6.1. Age-structured infection and disease development 
matrix. Here row values indicate infection age group (H(j)) and 
column values for age group in which infected individual devel- 
oped disease (D(j)). An individual who acquired the infection 
in j, and developed disease in j ; + u, is indicated by the cell 



more cells present where the condition j < f is satisfied, and also j is very 
low. (In fact, the average incubation period is not beyond a certain duration. 
It is not intended in the matrix to suggest that the lower the value of j then 
the larger the value of incubation period). If the age of infection is higher, for 
some j, and towards the last few possible age groups, then it is possible that 
f —j is shorter because individuals die naturally in old age. At the same time, 
the chance of infection in the very higher age groups (say 60+) is negligible for 
HIV (unless there are some rare causes). In the absense of age specific cohorts 
of infected individuals and follow-up data, it is not feasible to calculate disease 
progression rates and survival probabilities using direct cohort methods. In 
this section, we extend the method given in section 2 to estimate the average 
disease progression rates (or average incubation periods) for infections in age 
group j. This method is dependent on infection densities and data on disease 
occurances for the age group j. 

Let p(t,j) and q(t,j) be the probability density functions of infection density 
and incubation period for the age group j. If Q(t,j) is the distribution function 
of the incubation period, then Q(t,j) = q(t,j)dt. Now, the convolution of 
p(t,j) and Q(t,j) is given by 
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C(s,t) = / p{t,j)Q{t-s,j)ds. 

J —oo 

We call C the convolution of p and Q (i.e. p*Q, where * is the convolution 
operator). Therefore, 



P*Q = / p(t,j)Q(t - s,j)ds. 

J — oo 

Suppose an individual is diagnosed with a disease at age j in the year Uk- 
Then there is a possibility that this individual acquired the infection in any 
of the years prior to Uk (provided this individual is born in the year > Uq). 
Similarly, all those individuals who are diagnosed with the disease at age j + w 
in the year U n have actually acquired infection in any of the years from Uq to 
U n . In the same way, an individual infected at age j will be diagnosed with 
the disease in an age group > j. We consider model (16.11) . where four types of 
drugs were considered in section 2. 

Let Ao(j),Ai(j), A 2 (j) ,A 3 (j) be the parameter sets in age group j for the 
four kind of drugs. Let B (j),Bi(j), B 2 (j),B 3 (j) be the parameter sets C and 
E (j),Ei(j), E 2 (j),E 3 (j) be the corresponding events of diagnosis of disease 
in the age group j for the four types of drugs. The cumulative number of 
diagnosed disease cases up to U n for individuals who are diagnosed in the age 
group j is 



J(U < s < U n ,j) = J 0'*>j)> 

j*=o 

where I(0,j),I(l,j),...,I(j,j) are the numbers of disease cases diagnosed 
in age group j, and acquired the infection in the age group 0, 1, j. 

HO,]) = / p(t,0)Q(t-s,j)ds+ / p(t,0)Q(t-s,j)ds 
Jo Ju n . k 

[■Un-m l-Un 

+ p(t,0)Q(t-s,j)ds+ p(t,0)Q(t-s,j)ds 

J Un — l J Un — m 

p{t,0/A o )Q{t-s,j/B )ds+ / p{t,0/A 1 )Q{t-s,j/B 1 )ds 
+ p(t,0/A 2 )Q(t-sJ/B 2 )ds+ p(t,0/A 3 )Q(t-s,j/B 3 )ds 

J U n — I J Un — m 
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p U n _ k p U n — i 

1(1, j) = / p(t,l)Q(t-s,j)ds+ / p(t,l)Q(t-s,j)ds 
Jo Ju n _ k 

+ p(t,l)Q(t-s,j)ds+ p(t,l)Q(t-s,j)ds 

JU n -i JUn-m 
pU n — k pUn — l 

= / p(t,l/A )Q(t-s,j/B )ds+ / p(t,l/A 1 )Q(t-s,j/B 1 )ds 
Jo Ju n _ k 

fUn-m rU n 

+ p(t,l/A 2 )Q(t-s,j/B 2 )ds+ p(t,l/A 3 )Q(t-sJ/B 3 )ds 

JUn-l JU n - m 



= / p(t,j)Q(t-s,j)ds+ / p(t,j)Q(t-s,j)ds 
Jo Ju n - k 



I p(t,j)Q(t-s,j)ds+ / p(t,j)Q(t-s,j)ds 



rU n -k rU n -i 
= / p(t,j/A )Q(t-s,j/B )ds+ / p(t,j/A 1 )Q(t-s,j/B 1 )ds 
Jo Ju n _ k 

+ p(t,j/A 2 )Q(t-s,j/B 2 )ds+ p(t,j/A 3 )Q(t-s,j/B 3 )ds 

Similarly for unstructured populations, we assume that U n - k is the time 
of the introduction of drugs after the first year of detection of the disease in 
U . If E x {j) e [U n _ k , U n -{) and E 2 (j) e [U n _ u U„), then Z x {j) < Z 2 (j), 
otherwise if E 2 (j) e [U n _ k , U n -i) and E^j) E C/ n _ m ) then Z 2 (j) < 

Zi(j). If Exd), E 2 (j) e [U n - k , U n - m ), then Z^j) = Z 2 (j). Given an individual 
who was diagnosed with the disease in the age group j before U n is already 
developed in one of the four intervals [U ,U n - k ), [U n - k ,U n -i), [U n -i,U n - m ) 
and [U n - m , U n ). If E (j) e [L^/_i, U v ) C [U Q , U n _ k ), (for drug type i'), then the 
conditional probability of occurence of E (j) given E(j) is 
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Pr[E (j)/E(j)} = Pr[U i ,- 1 <J<U i ,,j/J<U n ] 



J [Aq(J), B (j)/Ui>,j] - J [Aq(J), B (j)/ Ul ,-uj] 
J[A (j),B (j)/U n ,j] 

where J values for E (j) are given by: 



J[Ao(j),B (j)/Ui,,j] = 



p(t,l/A )Q(t-s,j/B )ds 



/ p(t,0/Ao)Q(t-s,j/B o )ds 
Jo 

+ / p(t,j/A )Q(t-s,j/B )ds 
Jo 



J[A Q {j),B Q {j)/U i ,- U 3\ = [ p(t,0/A o )Q(t-s,j/B o )ds 

Jo 

/•fi'-i f u i'-i 
+ / p(t,l/A )Q(t-s,j/B )ds ■■■ +/ p{t,j/Ao)Q{t-s,j/B )d£ 
Jo Jo 



p{t,0/A )Q(t-s,j/B )ds 



+ / p(t,l/Ao)Q(t-s,j/B )ds ■■■ +/ p(tJ/A )Q(t-sJ/B )ds 
Jo Jo 

The above probability expressions are for the case without drug interven- 
tions. When drugs were initiated at U n -k, then these probabilities changed 
according to the occurence of Ei(j), E 2 (j), E 3 (j). Suppose Ei(j) D E 2 (j) = 0. 
If [E/fc'-i,!/*') C [U n - k ,U n -i), and E 1 £ [U n - k ,U n -i), then 



PriE^/EU)} = Pr[U kl -i<J<U kl ,j/J<U n ] 

J[A^ 3 %B^3)/Uy,j] -J[A 1 (j),B 1 (j)/u k/ ^j} 
J[A 1 (j),B 1 (j)/U n ,j] 



where J values for E\(j) are given by 
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J^UlB^yU^j] = / p&O^Qit-sJ/Bjds 

Jo 

+ f " p{t,l/A 1 )Q(t-s,j/B 1 )ds ■■■ + [ " p{t,j/A 1 )Q(t-s,j/B 1 )ds 
Jo Jo 

J[A 1 (j),B 1 (j)/U k ,- 1 ,j} = / p(t,0/A 1 )Q(t-s,j/B 1 )ds 

Jo 

fUy_ 1 f^k'-l 

+ / pfal/AJQit-sJ/Bjds ■■■ + p(t,j/A l )Q(t-s,j/B l )ds 
Jo Jo 

J[A l {j),B l {j)/U n ,j} = / pfaO/AJQit-sJ/Bjds 

Jo 



+ / p(t,l/A 1 )Q(t-s,j/B 1 )ds ■■■ + p(tJ/A 1 )Q(t-sJ/B 1 )ds 
Jo Jo 

In the above, instead of E 1 (j), if E 2 G [U n - k , U n ^i), then the probabilities 
would be 



Pr[E 2 (j)/E(j)] = Pr[U k ,^<J<U k ,,3/J<U n ] 

J[A 2 (j),B 2 (j)/U k ,,j] -J[A 2 (j),B 2 (j)/u k ,_ u j] 
J[A 2 (j),B 2 (j)/U n ,j] ' 

where J values for E 1 (j) are given by 



k> 



J[A 2 (j),B 2 (j)/U k ,,j] = / p(t,0/A 2 )Q(t-s,j/B 2 )ds 

Jo 

+ [ ' p{t,l/A 2 )Q(t-s,j/B 2 )ds ••• + [ ' p{t,j/A 2 )Q(t-s,j/B 2 )ds 
Jo Jo 

J[A 2 (j),B 2 (j)/U k ,- 1 ,j] = / p(t,0/A 2 )Q(t-s,j/B 2 )ds 

Jo 

+ / p(t,l/A 2 )Q(t-s,j/B 2 )ds ■■■ + p(t,j/A 2 )Q(t-s,j/B 2 )ds 
Jo Jo 
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J[A 2 (j),B 2 (j)/U n ,j} = / p(t,0/A 2 )Q(t-s,j/B 2 )ds 

Jo 

+ / p(t,l/A 2 )Q(t-s,j/B 2 )ds +/ p(t,j/A 2 )Q(t-s,j/B 2 )ds 
Jo Jo 

If [Ui-!,Ui) C [U n -i,U n - m ), and E x E [U n - h U n - m ), then 
Pr^OO/^O')] = Pr[U v ^<J<U v ,j/J<U n ] 

J[A 1 (j),B 1 (j)/U l ,,j] -J[A 1 (j),B 1 (j)/u l ^ u j] 
J[A 1 (j),B 1 (j)/U n ,j] 

where J values for Ei(j) are given by 

J[A 1 (j),B 1 (j)/U l ,,j] = f 1 ' p(t,0/A 1 )Q(t-s,j/B 1 )ds 

Jo 

+ f ' p(t,l/A 1 )Q(t-s,j/B 1 )ds ••• + / ' p(t,j/A 1 )Q(t-s,j/B 1 )ds 
Jo Jo 

J[A 1 (j),B 1 {j)/Ui l . 1 ,j] = / pfaO/AJQd-sJ/Bjds 

Jo 

+ [ l 1 p(t,l/A 1 )Q(t-s,j/B 1 )ds ••• + / ' ' p(t,j/A 1 )Q(t-s,j/B 1 )ds 
Jo Jo 

JlA^B^yU^j] = / pfaO/AJQd-sJ/Bjds 

Jo 

+ / p{t,l/A l )Q{t-s,j/B l )ds ••• + / ■p(t,j/A 1 )Q(t-s,j/B 1 )ds 
Jo Jo 

Suppose [Ui-!,Ui) C [[/„_/, C/ n _ m ), and £ 2 G [[/„_;, C/ n _ m ) then 

Pr^OO/^O')] = JV[t>i'-i< J<U v ,j/J<U n ] 

J [A 2 (j), B 2 {3)/U v ,j] - J [A 2 (j),B 2 (j)/ Ui ,. u j] 
J[A 2 (j),B 2 (j)/U n ,j] 
where J values for E 2 (j) are given as below: 
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J[A 2 (j),B 2 (j)/U l ,,j] 

+ / p(t,l/A 2 )Q(t-s,j/B 2 )ds 
Jo 

j[A 2 (j),fi 2 (j)M-i,j] 



+ 



Ui>. 



p(t,l/A 2 )Q(t-s,j/B 2 )ds 



J[A 2 (j),B 2 (j)/U n ,j] = 



/ p(t,0/A 2 )Q(t-s,j/B 2 )ds 
Jo 

+ / p(t,j/A 2 )Q(t-s,j/B 2 )ds 
Jo 

/ p(t,0/A 2 )Q(t-s,j/B 2 )ds 
Jo 

+ I 1 1 p(t,j/A 2 )Q(t-s,j/B 2 )ds 
Jo 

/ p(t,0/A 2 )Q(t-s,j/B 2 )ds 
Jo 



+ [ p{t,l/A 2 )Q(t-s,j/B 2 )ds ■■■ + [ p(t,j/A 2 )Q(t-s,j/B 2 )ds 
Jo Jo 

If E^j) = E 2 {j) e [£V-i,£/ P C [U n . k , U„) i.e. Z 1 (j) = Z 2 {j), then the 
conditional probabilities contain the same parameter sets. The probabilities 
for this situation are 



u„ 



Pr [E^j) = E 2 {j)/E{ 3 )] = Pr [U^ < J < U p ,,j/J < U n ] 

J[A 2 (j),B 2 (j)/U pl ,j] -J[A 2 (j),B 2 (j)/ Up/ ^,j] 
J[A 2 (j),B 2 (j)/U n ,j] 

where J values for E 2 (j) are given by 

J[A 2 (j),B 2 (j)/U p ,,j] = / p(t,0/A 2 )Q(t-sJ/B 2 )ds 

Jo 



+ / p(t,l/A 2 )Q(t-s,j/B 2 )ds 
Jo 

J[A 2 (j) : B 2 (j)/U pl ^j] 

+ / " 1 p(t,l/A 2 )Q(t-s,j/B 2 )ds 
Jo 



+ / p(t,j/A 2 )Q(t-s,j/B 2 )ds 
Jo 

/ p(t,0/A 2 )Q(t-s,j/B 2 )ds 
Jo 

+ / ' 1 p(t,j/A 2 )Q(t-s,j/B 2 )ds 
Jo 
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J[A 2 (j),B 2 (j)/U n ,j] = / p(t,0/A 2 )Q(t-sJ/B 2 )ds 

Jo 

+ / p(t,l/A 2 )Q(t-s,j/B 2 )ds ••• + / p(tJ/A 2 )Q(t-sJ/B 2 )dz 
Jo Jo 



Since Z 3 (j) > {Z (j), Z x (j), Z 2 (j)}, suppose E 3 (j) e U m >) C [U n - m ,U n ] 

now above probabilities are 



Pr[E 3 (j)/E(j)} = Pr[U m .- X <J <U m ,,j/J <U n \ 

J[A 3 (j),B 3 (j)/Um>,j] -J[A 3 (j),B 3 (j)/u m ,^j] 
J[A 3 (j),B 3 (j)/U n ,j] 

where J values for E 3 (j) are given by 

ru m > 

J[A 3 (j),B 3 (j)/U m ,,j] = / p(t,0/A 3 )Q(t-s,j/B 3 )ds 

Jo 

r U m> f U m> 

+ / p(t,l/A 3 )Q(t-s,j/B 3 )ds ■■■ + p{t,j/A 3 )Q{t-s,j/B 3 )ds 
Jo Jo 



J[A 3 {j),B 3 {j)/U m ,- 1 ,j\ = [ m ~ 1 p(t,0/A 3 )Q(t-s,j/B 3 )ds 

Jo 

f U m'-l f U m'-l 

+ / p(t,l/A 3 )Q(t-s,j/B 3 )ds ••• + / p(t,j/A 3 )Q(t-s,j/B 3 )ds 
Jo Jo 

J[A 3 (j),B 3 (j)/U n J] = / p(t,0/A 3 )Q(t-sJ/B 3 )ds 

Jo 

+ / p(t,l/A 3 )Q(t-sJ/B 3 )ds ••• + / p(tJ/A 3 )Q(t-sJ/B 3 )ds 
Jo Jo 



Using the above conditional probabilities, likelihood functions are constructed 
by assuming some parametric form for the diagnosed disease cases. For each 
age group above, analysis is conducted to estimate the incubation periods by 
age group. 
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7. Conclusions 

The improved models that address the impacts of anti-retroviral therapy, 
protease inhibitors and combination of drugs presented in section 1 seem useful 
in understanding the dynamics of variables for individuals with the full blown 
disease for no-drug, drug 1, drug 2 and drug 3, i.e D ZQ , D Zl , D Z2 and D Z3 . 
Using the methodology in sections 2 to 4, (despite being lengthy), we are 
able to estimate the parameters for the incubation period for each drug type, 
by the deconvolution method. We have demonstrated this method for three 
types of drugs, and one can obtain B for as many drugs as possible from the 
formulae for n— types of drugs in section 3. So far, there is no evendence of 
drugs being useful in avoiding contracting the disease. Drugs may be useful for 
avoiding opportunistic infections for some specific periods of time. Eventually, 
an individual will succumb to AIDS, whether or not that individual takes drugs 
(which is also demonstrated in the truncation effect). The truncation effect 
formulae can be used to obtain the parameter set (say, B T ), but we did not 
demonstrate this numerically. 

We did not introduce intracellular delay that might arise due drug interven- 
tions. There are not many quantititive results available on the relationship 
between the dose of a drug and the resultant delay in the development of 
the disease. Suppose s±, S2, S3..., Sk are k levels of doses of a single drug, and 
Ti, T2,T3, ...,Tfc are the respective delays obtained in producing a new infected 
cell. Then we can write the relation R 2 (s, r) between s and r as 

{ELi fa-?)} 

{ELi(^-*)} 2 {E?=i(^-r)} 2 

R (s, r) is called the correlation coefficient of dose-delay, s is the mean 
dose-level and r is the mean delay. This experiment can be conducted for 
various doses (say) for drug type j = l,2,3...n. Each drug will produce 
a delay depending upon the dose level. From this, the average delay can be 
statistically compared to understand the mean dose effect due to a particular 
drug, and hence the drug efficacy. However, this does not give dynamics over 
the time period, but it is very useful in preparing the baseline parameters for 
simulation studies, and also for the models explained in sections 1, 2 and 5. 
There might be a possibility of exploring the impact of delay in the conditional 
probabilities expressed in this work. 

Our work may be interesting for people working on developing computa- 
tional techniques for solving integro-differential equations, algorithms to solve 
convolution type equations in epidemiology, and EM-type algorithms. The 



RUNNING HEAD: VARYING INCUBATION MODEL 



35 



age-structure analysis presented is more complicated than analysis presented 
for the non-age structured populations, and we provide a new kind of analysis 
for the incubation period. When reported disease cases and densities of the 
infection are available for a period of several years in the population, then 
this kind of analysis offers a reliable method to estimate the incubation period 
distribution. 
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Appendix I : Conditional probabilities for generalized 

MULTIPLE DRUG IMPACT 

Here we derive expressions for conditional probabilities when several drugs 
are available, and the incubation period is non-monotonic. When such a sit- 
uation arises there will be several combinations of orders of Zs. We take one 
such situation and write corresponding Ls for the purpose of demonstration. 

Suppose Z < ... < Zk = ... = Zk +n+ i < ... < Z^j. Let us divide this into the 
following two inequalities and an equality: Z$ < ... < Z^ , Zk+i = ... = Zk +n 
and Zk+n+i < ■■■ < Zn- If we consider the first and third inequalities, then 



u Nq 



V(A,B/U Nk ) = / h(t/A No )G(t-s/B No )ds + 
Jo 

fU Nl 

/ h(t/A Nl )G(t- s/B Nl )ds 
Ju No 

r UN k 

•••+ / h(t/A Nk )G{t-s/B Nk )d& 



U N n+k+1 



V(A,B/U Nn ) = I h(t/A Nn+k+1 )G(t-s/B Nn+k+1 )ds + 

o 



f UN n + k + 2 

/ Ht/A Nn 

+ k+2 

)G(t - s/B Nn+k+2 )ds 

Ju No 

r u ^ N 

■■■ + h(t/A NN )G(t- s/B NN )ds 



We can express {P (E Ne / E)} d g ~Q and {P (E Nf) / E)} 9 9 ~^ +k+1 and the correspond- 
ing {L^ e } g as shown in the section 3. Then L^ e is maximised for the set 
[Ag,^]. We obtain N — n — k sets of [A, B] values, and the corresponding 
likelihood functions {L Ne } g=0 and {L Ne } g=n+k+1 
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Figure 7.1. Truncated incubation period. The idea of trun- 
cated cumulative distribution of the incubation period is plot- 
ted. After a certain time duration, there will not be any gain 
due to therapy. Median incubation period is represented by the 
line cutting the curve at 0.5, corresponding to the Y-axis. 



Appendix II: Truncated incubation period 

Suppose there is an upper bound for the impact of drugs on the incubation 
period; that is, the incubation period cannot be increased after a certain time 
point after the drug use. Then the likelihood equations explained in section 
4 would change accordingly. There was an attempt earlier to truncate the 
incubation period with the help of the truncated Weibull distribution (Rao 
and Kakehashi, 2005). They have not seen the impact of drugs using such 
functions. If Z, the length of the incubation period, and if Z c is the truncation 

point, then G(Z) = 1 - exp j- ^)* 2 J , for < Z < Z c , and G(Z) = 

SlJ j exp j- (|) (|) (fe ~ 1)(2 " c) j, for Z > Z c . Here, S 1 ,5 2 

are scale and shape parameters. One can construct a likelihood function for 
each drug type using such functions as follows: 



1 - exp <; - i 4-\ 2 



(7.1) 



L(A,BAR,) = L <Zc + L> Zc 



where 
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l< Zc = YihUfaw-Yihu-vw) 



and 



biU) 

h(j - 1) 

b\(j) 



b\(j ~ 1) 



bl(j) 



(is 



e Qls2+a2S+Q3 ^ 1 - exp I - ( |- 



<5a' 



ds 



eaiS 2 +Q2S+as I 1 _ ^ 



9 I / 2 



! 1 — exp < — 



Si 



exp 



SA (tc 

Si U 



S-2 ' 



(<5 2 -l)(z-z c )' 



ds 



1 — exp 



x 



exp < - 



S 2 \ ft c 



Si \S 



(s 2 -i){z- Zc y 



ds 



1 — exp < — 



Si 



&2 



x 



exp 



SA (U 
Si) Ui 



(<5 2 -l)(z-z c )' 



ds 



For each drug type, an expression of the type 17.11 can be derived. Despite 
the assumption on truncation as mentioned above, the incubation period could 
vary according to the type of the drug. 



